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enables the design of switchable spin-based quantum gates by addressing only the electron spin. 
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Since the idea of quantum information processing 
(QIP) fascinated the scientific community []], [3,0] 5 elec- 
tron and nuclear spins have been regarded as promising 
candidates for quantum bits (qubits) [4]. Nuclear mag- 
netic resonance (NMR) in liquid state was the first spec- 
troscopic technique used to demonstrate several quantum 
computation algorithms 

H S0, i, while the emerged 
scalability limitations can be overcome by using_electron 
paramagnetic resonance (EPR) spectroscopy [9j]. The 
construction of quantum gates based exclusively on elec- 
tron spins with controllable exchange interactions is how- 
ever a challenging task, especially when crystalline ma- 
terials of this kind are on demand. This difficulty has 
motivated an on-going effort to find appropriate elec- 
tron spin qubits, including for instance quantum dots 
[Tot , single- molecule mag nets fll| , and antiferromagnetic 
heterometallic rings [12[ or metal clusters [13]. On the 
other hand, hybrid electron-nuclear spin systems with 
long decoherence times can be found in a variety of mate- 
rials like organic single crystals [3] , endohedral fullerenes 
Til . T3 , 17]. phosphorous donors in silicon crystals [lil ]. 
and nitrogen- vacancy centres in diamond [19, 20]. They 
have been used to perform two-qubit quantum opera- 
tions, demonstrate entangled states, or build solid-state 
quantum memories. These systems benefit from well- 
defined and separated EPR and NMR transitions that 
can be selectively manipulated by resonant mw and ra- 
dio frequency (rf) radiation, respectively. Therefore, ad- 
vanced EPR methods [21] employing both selective mw 
and rf pulses play a key role in QIP based on electron- 
nuclear spin systems. 

The significant difference between decoherence times of 
electron (T 2e ) and nuclear (T 2n ) spins, which is due to 
the fact that their gyromagnetic ratios ji differ by two to 
three orders of magnitude, has been utilized for the real- 
ization of a quantum memory that uses slow relaxing nu- 
clear spins to store the information and fast relaxing elec- 
tron spins for processing and readout [18[. On the other 
hand, this difference also restricts the clock rates of gate 
operation to the corresponding ones of the slow relaxing 
qubit, namely to the MHz frequency scale. More seri- 



ously though, the combination of the small gyromagnetic 
ratios 7 n of nuclei with the currently available rf fields B\ , 
results in slow rotations of nuclear spins. For instance, 
the nutation (Rabi) frequency of a proton nuclear spin for 
a typical value of B\—\ mT is uj\/2it = j n Bi/h =42.6 
kHz which implies a length of At^ = 11.7/is for an rf 
7r-pulse. This time interval is of the same order of mag- 
nitude with typical electron decoherence times, T^e, and 
thus the efficiency of the two-qubit quantum gate may 
become questionable due to relaxation losses. 

One way to overcome this difficulty is to utilize the 
hyperfine interaction which is typically stronger than the 
Rabi nuclear frequency, i.e. of the order of some MHz 
in organic radicals. When this interaction is anisotropic 
the quantization axes of the hyperfine fields deviate from 
the z axis (see FIG. QJi) and thus the nutating nuclear 
spin can be inverted by implementing a proper set of 
mw 7r-pulses separated by suitable time delays Q • This 
concept has been recently used in a similar fashion [2^ 




FIG. 1: (Color online) (a) Bloch sphere showing the trajec- 
tories of nuclear spin magnetization under free evolution as- 
suming that the nuclear spin is polarized in the beginning 
(jm=-\- 1/2). The quantization axis of the effective magnetic 
field experienced by the nucleus depends on the state of the 
electron spin: for ms=+l/2 the nuclear spin precesses about 
hex (yellow trace) with angular frequency oo a = |o^i2 1 ; for 
ms=-l/2 the nuclear spin precesses about ftp (green trace) 
with angular frequency up = \uj3a\- (b) Exact cancellation 
case, uji — A 1 '2. 
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FIG. 2: (a) Schematic representation of the stable radical 
•CH(COOH)2 and the principal axes system of the hyperfine 
tensor, (b) Energy level diagram depicting the allowed (al, 
a2) and forbidden (fl, £2) electron spin transitions, (c) EPR 
stick spectrum. 



to demonstrate the creation of nuclear spin coherence via 
amplitude- modulated mw pulses [24]. 

Herein we focus on a special case of electron spin 
echo envelope modulation (ESEEM) spectroscopy, the 
so-called exact cancellation [2l| where the hyperfine and 
electron Zeeman interaction cancel each other in one ms 
manifold. FIG. [TJd shows that in this case the effective 
magnetic field at the nucleus is either almost parallel or 
exactly perpendicular to the z axis, depending on the 
state of the electron spin. Consequently, during free 
evolution of the system with ms = —1/2, the nuclear 
spin naturally oscillates between the raj = +1/2 and 
raj = —1/2 states without the need of rf field. In ad- 
dition, it can be locked in a desired state by inverting 
the electron spin state to ras=+l/2 with a mw 7r-pulse 
applied after a proper time delay. 

To demonstrate this control, a single crystal of 7- 
irradiated malonic acid is used. Here the paramagnetic 
species is the stable radical CH(COOH)2 shown in FIG. 

where the unpaired electron (5=1/2) resides on the 
carbon 2p z orbital and is hyperfine-coupled to the a- 
proton nuclear spin (i = l/2)[25j. The rotating frame spin 
Hamiltonian is given by [2l| 



Ho = V S S Z + u T I z + AS Z I Z + BS Z I X 



(1) 



where Qs = — u mw is the offset of the electron 
Zeeman frequency uos = g/3 e Bo/h from the mw fre- 
quency uOrnwi wi = — g n /3 n Bo/h is the nuclear Zeeman fre- 
quency, g and g n are the electron and nuclear g- factors, 



(3 e and (3 n are the Bohr and nuclear magnetons, B Q is 
the static magnetic field along £-axis, and A, B describe 
the secular and pseudo-secular part of the hyperfine cou- 
pling. In this four-level electron-nuclear spin system 
(FIG. [2}d) there are six possible transitions: four EPR 
with Aras = ±1 as illustrated in the stick spectrum of 
FIG. [2t, and two NMR with frequencies 0012 and 0034. 
For a proton nuclear spin (cjj < 0) and a negative hy- 
perfine coupling (A < 0) in the weak-coupling regime 
(\A\ < 2|c^j|), the singed nuclear transition frequencies 
can be expressed as 0012 = — \J (u) 1 + A/2) 2 + (B/2) 2 
and 0034 = -vVj - A/2) 2 + (B/2) 2 . In addition, the 
quantization axes of the effective frequency vectors are 
determined by the angles r] a = arctan[— B/(A + 2cjj)] 
and rip — arctan[— B / (A — 2a; j)] which are defined rel- 
ative to z direction. For the case of interest of exact 
cancellation, A = 2o;j, the nuclear transition frequen- 
cies become LO12 = — y (2c<;j) 2 + (B/2) 2 and CJ34 = —B/2 
(with B > 0) with the corresponding angles being rj a = 
arctan[— B/Auoi) and rjp = —tt/2. 

The concept of manipulating the proton nuclear spin 
solely by using mw pulses can be best described in terms 
of the density matrix formalism. Figure [3K shows an ex- 
ample of this control assuming that the system is initially 
prepared in the pseudopure state <jp a = \/3a)(/3a\ which 
is described by the density matrix in the product basis 
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During free evolution of the system under the spin Hamil- 
tonian (Eq. [T]) the hyperfine field drives the nuclear spin 
transition and the population is periodically transferred 
from ap a to a /3 /3 and vice- versa with frequency CJ34. This 
transfer is complete (i.e. a (3/3 = \(3f3){f3l3\) only if the 
exact cancellation condition sin rjp = — 1 is fulfilled, and 
occurs for free evolution times t\ — (2ra + /uj^ (ra 
( ) . 1.2....) (see supplementary information). A subse- 
quent semi-selective mw 7r-pulse, that simultaneously ex- 
cites the forbidden 23 and allowed 24 EPR transitions 
[H, (represented as (7r)2324), transfers the popula- 
tion to the I a(3) state. Provided that the angle r/ a is 
very small, the populations and nuclear coherences in 
the ras = +1/2 subspace are virtually trapped because 
of the negligible branching between \aa) and \a(3) states. 
Under that treatment the nuclear spin can be locked in 
the raj = —1/2 state for arbitrary time, T2. An addi- 
tional (7r)2324-pulse will transfer back the population to 
the \/3/3) state and the nuclear spin will start oscillating 
again between raj = +1/2 and raj = —1/2 with fre- 
quency CJ34. Quantitatively, this is demonstrated in FIG. 
[3)3 that shows the simulation of the nuclear spin polar- 
ization, (I z ). Therefore, by choosing appropriate time 
delays and semi-selective mw pulses, the full control of 
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FIG. 3: (a) Representation of the density matrix in the product basis |aa), \ctf3), \/3a), \/3/3) during nutation or locking time 
periods, (b) Numerical simulation of the nuclear spin polarization (I z ) during the applied pulse sequence that uses semi-selective 
mw 7r-pulses (shown at the top). Spin Hamiltonian parameters: nuclear Zeeman frequency, uoi/2n = —14.58 MHz; hyperfine 
coupling constants, A/2tt — —29.06 MHz, B/2tt — 6.45 MHz. (c) Corresponding numerical simulation of the electron spin 
coherence (S x ) after the detection pulse sequence (71-/2)2324 — r — (7r)2324 — t— echo with r = 2ti /uosa =310 ns. 



the nuclear spin can be achieved. 

Although direct measurement of nuclear spin polariza- 
tion is not possible in EPR spectroscopy, FIG. shows 
that the electron spin coherence (S x ) after the sequence 
(7r/2) 2 324-T-(7r)2324-7--e-cho with r 27rra/cj 34 , (ra 
1,2,3...) represents the population difference between 
states 2 and 4. Consequently, such a sub-sequence can 
be used to monitor the state of our system. 

For the experimental demonstration of this control the 
system can start from a pseudopure state, e.g. crpa- For 
the case of isotropic hyperfine coupling (5 = 0), such a 
state can be prepared with a selective mw pulse P24(A))? 
with /?o = arctan(-l/3) = 109.5°, followed by a tt/2- 
rf pulse, Pi2(7r/2)|4|. For anisotropic hyperfine coupling 
and a strong mixing of states the above pulse sequence 
does not result in the desired pseudopure state. However, 
for the special case of exact cancellation, the operation of 
a selective mw pulse P2a(P) (with effective rotation angle 
(3 e ff = arctan(— 1/3) = 109.5°) at the thermal equilib- 
rium density matrix, ctq = — S z , gives the pseudopure 
state a aa in good approximation (see supplementary in- 
formation). An additional semi-selective (tt)i314 pulse 
applied after time > 5T2 e complements the prepara- 
tion part of the pulse scheme shown in FIG. [4^ and leads 
the system to the pseudopure state a pa. 

The experiments were performed at a crystal orienta- 



tion for which the hyperfine coupling is very close to ex- 
act cancellation: the three-pulse ESEEM spectrum (FIG. 
Ht) measured at the observer position Bo = 343.7 mT 
(uji/27t = —14.6 MHz) shows two peaks at frequencies 
cji 2 /2tt = -29.42 MHz and ^ 34 /27r = -1.83 MHz from 
which the angles r] a = —3.6° and rjp = —87.1° can be in- 
ferred. Figure [4)3 shows that the two EPR transitions 23 
and 24 are not resolved due to the small nuclear frequency 
and inhomogeneous broadening. Therefore, an accu- 
rate excitation of transition 24 with a selective P24 pulse 
is difficult to be accomplished. Moreover, the necessity 
of the semi-selective (71") 1314 pulse requires an additional 
mw frequency which was not available with the current 
experimental setup. However, both technical issues were 
circumvented by replacing the preparation part of FIG. 
2k with a semi-selective (7r)2324 pulse that interchanges 
populations between levels 2 and 4. Their population 
difference (which is our observable here) is expected to 
oscillate between zero and a maximum value because af- 
ter the semi-selective pulse levels 3 and 4 will start ex- 
changing populations, in analogy to the situation of pop- 
ulation transfer during free evolution of crpa under the 
spin Hamiltonian (Eq. [Tj). Figure HJi (gray trace) shows 
that the signal virtually vanishes for t = (m + 1/2) T 
(m = 0, 1, 2, ...), where T = 504 ns (« 271-/^34), implying 
that our system is indeed very close to exact cancella- 
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FIG. 4: (a) The proposed pulse sequence consists of three 
parts: preparation, for creating the pseudopure state ap a ] ma- 
nipulation, consisting of a time delay r n and a semi-selective 
(tt)2324 or (7r)i3i4 pulse; detection of electron spin coherence 
through a two-pulse echo with r = 2tv/oo34- (b) Free induction 
decay (FID)-detected EPR spectrum. The arrow indicates the 
observer position Bo = 343.7 mT of the semi-selective (7r)2324 
pulses, (c) Three-pulse ESEEM spectrum, (d) Electron spin 
echo (ESE) intensity as a function of time tdetection between 
preparation and detection. 



tion. By applying semi-selective mw pulses, the nuclear 
spin evolution can be locked in the mj = — 1/2 state or 
released at will (FIG.[Ui, black trace). 

In conclusion, we have demonstrated that our approach 
provides a fast way of controlling nuclear spins in hybrid 
spin systems and thus overcomes the problem of relax- 
ation incompatibility between electron and nuclear spins. 
Furthermore, the method offers a means for controlling 
the interaction between the two qubits in a fast and ef- 
ficient way. Although the hyperfine interaction is always 
present, the special case of exact cancellation considered 
here induces two different effective states of interaction 
regarding the behaviour of the nuclear spin: an active 
(nutation) and a passive state (lock) that can be switched 
with a single mw pulse. While the control of dipolar in- 
teractions between neighbouring electron spin qubits has 



been proposed using SWAP operations [28|, |29[, to the 
best of our knowledge to date this kind of control in sim- 
ilar spin systems has been experimentally demonstrated 
only once using decoupling methods [30[. 

Our work shows that advanced pulsed EPR methods 
like ESEEM spectroscopy can play a key role in QIP 
based on hybrid electron- nuclear spin systems. Further- 
more, it gives new perspective to such systems by consid- 
ering them not only for performing quantum memories 
but also for building solid-state quantum gates. Given 
that the development of new materials as well as new 
methods is equally important for the realization of a 
quantum computer, it is likely that the scheme presented 
here will become a standard method for manipulating 
electron and nuclear spin qubits and controlling their hy- 
perfine interaction in the future real quantum computers. 
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THEORY 

A. Population transfer between \j3a) and \j3j3) as driven by the hyperfine field. 



The pseudopure state ap a is described by the density matrix in the Cartesian product 
basis (CPB) 

/o o o o\ 



°(3o 



(i) 




10 
\0 0/ 

The spin Hamiltonian Hq = Q,sS z + ujiI z + AS Z I Z + BS Z I X is diagonalized {Hq = UTLoU^) 
by the unitary transformation 

\ 



U = 



( cos(^) -sin(^) 
sin(^) cos(^) 



cos(f) -sin(f) 







(2) 



sin(f ) cos(f) / 

where rj a = arctan[— B/(A + 2a;/)], rj^ = arctan[— B/(A — 2a;/)], and rj = (rj a — r}p)/2. 
Therefore, the initial density matrix o\ = crg Q is expressed in the eigenbasis as 
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which, after free evolution under the spin Hamiltonian 7Yq for time t, becomes 



a' t = exp(-iHot)a[ exp(i^t) 
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(4) 



where 6J34 = —\J (uj — A/2) 2 + (B/2) 2 is the signed nuclear transition frequency in the 
^-electron spin manifold. This density matrix in CPB reads 
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^sin(2^) sin 2 (^) + 
- sin(^) sin^t) 



\ sin 2 {rip) [1 - cos(w 34 t)] 



/ 



(5) 



From Eq. [5] we see that populations and nuclear coherences are periodic functions of time 
with angular frequency CJ34. The initial state a\ = ap a recurs for times T = 2m7r/uJs4 
(m = 1,2,3,...). At exact cancellation (sin(^) = —1 ) and times r = (2m + 1)tt/ujs4 
(m = 1, 2, 3, ...) complete population transfer occurs, i.e. a T = app. 



B. Preparation of the pseudopure state ap a . 

Starting from the thermal equilibrium density matrix ctq = — S Z) the standard method 
for preparing the pseudopure state ap a is the application of a selective mw pulse ^(A)), 
with nominal rotation angle f3 = arctan(— 1/3) = 109.5°, followed by a selective rf pulse, 
Pi 2 (7r/2). This is true only for isotropic hyperfine coupling, where 77 = 0. Here we examine 
the case of a non-zero anisotropic hyperfine coupling and particularly the case of exact 
cancellation. A selective mw pulse along the rotating frame y-axis, which is on-resonance 
with the allowed EPR transition 24 and has a nominal rotation angle /?, is given by 

= exp[-i/3 cosrjSy 4 ], with 5^ 4 = ^S y - S y I z (6) 



2 



Therefore, after this pulse the density matrix is expressed in the eigenbasis as 

\ 

-| cos 20 -± sin 20 

± 

\ -\ sin 20 \ cos 20 ) 

with 20 = (5 cosr}. When the condition of exact cancellation is fulfilled (sin^) 
transformation back to CPB gives 
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-1 ), the 
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4 cos 20 
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= 1 + ^08 20^/ 

(8) 



For clarity, the off-diagonal elements of the density matrix in Eq. [8] have been omitted. 
They represent electron and nuclear spin coherences that have completely decayed after 
time Td > 5T 2e . We can now determine the angle = | cost] in order to fulfil the condition 
B=C=D, i.e. to get the pseudopure state a aa . From this requirement we find 

^ , - COS 20 cos — - 
2 7 2 r v 2 ; 4 



sin ( — ) 



cos 20 



which gives 



cos 20 



-1 + 2sin 2 (^) 
2cos 2 (f ) + l 



1 
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6 



(9) 



The approximation in Eq. [9] holds for small angles rj a = arctan(— B/Aluj)^ i.e. when B <C 
4a; j. Therefore, application of a selective mw pulse P^iP) with effective rotation angle 
f3 e ff = /3cos7y = arctan(— 1/3) = 109.5° at the thermal equilibrium density matrix, gives in 
good approximation the pseudopure state a aa without the need of rf pulse. An additional 
semi-selective (7r)i3i4 pulse applied after time drives the system to the pseudopure state 
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II. MATERIALS AND METHODS 



Single crystals of malonic acid were grown from saturated aqueous solution by slow evap- 
oration and then irradiated at room temperature with a 60 Co 7-ray source. The occurring 
paramagnetic species is the stable radical -CH(COOH) 2 shown in Fig. 2a, where the un- 
paired electron (S = 1/2) resides on the carbon 2p z orbital and is hyperfine-coupled to the 
a-proton nuclear spin (/ = 1/2). The hyperfine interaction is highly anisotropic and exhibits 
rhombic symmetry with principal values {A X) A y) A z ) = (—29, —91, —61) MHz. 

Pulsed EPR measurements were performed at room temperature using a Bruker ESP380E 
spectrometer operating at X-band frequency (u mw /27r = 9.646 GHz). The field-swept EPR 
spectrum of Fig. 4b was recorded by integrating the FID after a selective mw tt/2 pulse 
of 500 ns. The three-pulse ESEEM spectrum was obtained with the pulse sequence tt/2 — 
r — tt/2 — T — tt/2— echo, with a tt/2 pulse of length 16 ns and an interpulse delay r = 200 
ns. Time T was increased in steps of dt = 8 ns starting from T = 56 ns (512 points). A 
four-step phase cycling was used to remove unwanted echoes. The time trace was baseline 
corrected with a two-order exponential, apodized with a Gaussian window, and zero-filled. 
After Fourier transform the absolute- value spectrum was calculated (Fig. 4c). 

Our experiments for the demonstration of the control of the proton nuclear spin were 
performed on the low field line of the EPR doublet as indicated by the arrow of Fig. 4b. 
7r-pulses of length 32 ns were used corresponding to mw field strength of Ui/27T = 15.6 
MHz. This is smaller than the hyperfine splitting at this crystal orientation (^29 MHz) and 
therefore the pulse excites only transitions 23 and 24 (which are separated by Us4/2tt = —1.83 
MHz). Such a pulse is called semi-selective and is represented as (7r) 2 324- The state of 
the system was monitored by measuring population differences with the pulse sequence 
(7172)2324— t— (7r) 2 324— t— echo that uses r = 552 ns. A four-step phase cycling was applied to 
remove unwanted echoes created by the semi-selective (7r) 2 324-pulses of the manipulation part 
of the control scheme. In Fig. 4d the preparation part of the pulse sequence is replaced by a 
semi-selective (7r) 2 324 pulse. After this pulse the signal oscillates with frequency ujp /2tt = 1.83 
MHz under free evolution (gray trace). By applying semi-selective (7r) 2 324 pulses at positions 
indicated by arrows the nuclear spin can be locked in the mi = —1/2 state (black trace, 
locking time intervals r 2 = 2016 ns and T4 = 1288 ns) and then released at any point (black 
trace, nutation time intervals r 3 = 1504 ns and r 5 = 2096 ns). 
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